Angular Momenta and Spin-Orbit Interaction of Nonparaxial Light in Free Space 
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I. INTRODUCTION 

The problem of the identification of the spin and or- 
bital parts of the angular momentum (AM) of an electro- 
magnetic wave has a long history and has posed funda- 
mental difficulties in both quantum electrodynamics and 
classical optics [H-ll]- 

It is knovifn that the photon AM operator in the mo- 
mentum (plane- wave) representation has the form [l[ : 



(k X ^k) + S = L + S . 



(1) 



Here the orbital part is L = f x p (p = k, f = i9k, 
k is the wave vector, and we use units h — c ^ 1), 
whereas S is the spin-1 operator given by 3 x 3 matrices 



(5a) 



{caij is the Levi-Civita symbol) that act 



on the Cartesian components of the wave electric field. 
Canonical orbital AM (0AM) and spin AM (SAM) op- 
erators, L and S, satisfy so(3) algebra and generate ro- 
tations in spatial and polarization degrees of freedom, 
respectively. However, "t/ie separation of the total AM 
into orbital and spin parts has restricted physical mean- 
ing. ... States with definite values of 0AM and SAM do 
not satisfy the condition of transversality in the general 
case.^^ In 1994, Van Enk and Nienhuis put forward 
an alternative, non-canonical AM separation, where the 
modified spin and orbital parts are measurable and con- 
sistent with the transversality of the wave, although they 
are not generators of rotations 

In classical optics, the two parts of Eq. ([T]) can be 
unambiguously associated with the QAM and SAM for 
paraxial light, where the eigenmodes of Lz = —id^ {4> 
is the azimuthal angle in k space) and Sz are circularly 
polarized vortex beams with the corresponding quantum 
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numbers i = 0,±1,±2, ... (topological charge of the vor- 
tex e'^"^) and a = ±1 (helicity) However, for non- 
paraxial fields the identification of 0AM and SAM meets 
serious difficulties 0-13 • Calculations based on the re- 
cently suggested division of the Poynting energy flow into 
spin and orbital parts show that the non-paraxial 

correction to the 0AM is proportional to cr rather than to 
^ 0, 01- This resulted in the conclusion that "m the gen- 
eral non-paraxial case there is no simple separation into 
£-dependent orbital and a-dependent spin component of 
AM" 0. 

In this paper we re-examine the problem and give an 
exact self-consistent solution in terms of both the fun- 
damental photon operators and classical energy flows. 
The identification of the well-defined measurable 0AM 
and SAM of light is shown to be closely related to the 
analogous problem for the position of localized photons 
[gt-fTTj. Our approach generalizes and unifies previously 
disjointed results: (i) non-canonical 0AM and SAM op- 
erators obtained earlier for the second-quantized fields 
0; (ii) non-commutative photon position operator and 
Berry monopole field in momentum space [lo| : and 
(iii) separation of the spin and orbital parts of Poynt- 
ing energy flows 0-0]. We find that the cr-dependent 
non-paraxial part of the 0AM arises from Berry-phase 
terms describing the spin-orbit interaction (SOI) of light. 
A similar effect occurs dynamically upon spin-to-orbital 
AM conversion in focusing and scattering of polarized 
light fq. [l^jlii . Other manifestations of the SOI are the 
spin |15l4l8l | and orbital [19h22l | Hall effects of light (i.e., 
£- and cr-dependent transverse shifts of the field center 
of gravity) that are described by ot ir p osition operator 
and take place even in free space [3 We apply 

the general theory to vector Bessel beams, for which the 
fundamental operators manifest themselves in immedi- 
ately observable i- and a-dependent intensity distribu- 
tions. The exact wave results are also explained in terms 
of the underlying geometrical-optics rays and caustics. 
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II. OPERATOR FORMALISM 

We consider an electromagnetic field in free space, 
characterized by its plane-wave electric-field spectrum 
E (k) without evanescent modes. The SOI of light origi- 
nates from the transversality constraint, k - E ~ 0, which 
couples polarization to the wave vector and reduces the 
full 3D vector space of the electric field components to the 
2D subspace of the components tangential to a sphere of 
directions in k-space. The operators L and S do not keep 
this subspace invariant, i.e., their action on a transverse 
mode results in a non-zero longitudinal component [U, . 
However, this subspace is invariant for the total AM op- 
erator J, and one can divide it into two parts consistent 
with the transversality condition: 

J = L'-^S', L' = t-KX (^KxSj , S' = • s) , (2) 

where k = k/fc and the modified 0AM and SAM op- 
erators L' and S' can be regarded as projections of the 
operators L and S onto the transversality subspace [1]. 

The modified SAM operator S' is proportional to the 
helicity operator a — k ■ S, whereas the 0AM operator 
can be written as L' = f' x k with 



k X S ^ k X S 



fc2 



(3) 



The modified position operator ([3]) has been considered 
in the context of photon localization and Berry phase 0- 
[T]| . It describes the observable center of gravity of the 
field and brings about the space non-commutativity with 
the monopole term in k-space: 



fc3 



(4) 



The operators L' and S' do not satisfy the so(3) AM 
algebra and have unusual commutation relations: 



= 0, 



^11 L'l 



Tl ql 



(5) 

At the same time, the modified operators transform as 



vectors under rotations; 



J.,0' 



o;, o' = L', s'. 



and f'. The commutation relations (O unveil the sim- 
ilarity of operators L' and S' to those obtained for the 
second-quantized fields in Q . Although they do not gen- 
erate rotations, it is suggested that they do correspond 
to observable continuous values of the 0AM and SAM of 
a non-paraxial transverse field • 

Remarkably, in the helicity representation the matrix 
components of the operators ^ and ^ become diag- 
onal. We introduce spherical coordinates {9, cj), k) with 
basic vectors {ee,e^, k) in k-space, so that the free elec- 
tric field has only {eg, e^)-components. The helicity basis 
of circular polarizations corresponds to the basic vectors 
6=*= = e='=""'^ {eg ± ie^) /a/2, where 6=*=*""^ is an arbitrary 



gauge factor [ll|. Transition from the global Cartesian 
field components {Ex,Ey, Ez)^ to the helicity amplitudes 
)'^ is realized via the local unitary trans- 
formation U {( 



Ra {a) 



<p) = Rz {-(f)) Ry {-9) {m<j)) V, where 
is the matrix of rotation by an angle a 
with respect to the a-axis, whereas V is the constant 
transformation from linear- to circular-polarization ba- 
sis. Making the transformation of operators ^ and ([3]) 
to the hehcity basis, 6' U^O'U, we obtain: 

S' ^ K.a , L' = —ik X 9k — As X k , (6) 
f' = idu — A.B , p = k , w = oj (7) 

Here we included the momentum and energy opera- 
tors, p and w (which are unaffected by the transfor- 
mations), uj is the frequency, the helicity is diagonal: 
(T = diag (1,-1,0), and 



k X S 



cos 9 



ksm9 



-ae^ 



(8) 



is the Berry gauge field (connection) which corresponds 
to the monopole curvature Fb = dk x Ab = ak/k^ 
[lol [TTj |. Hereafter we choose the gauge m = 1, which 
corresponds to the absence of the phase singularity (Dirac 
string) along the positive z-axis in Eq. (jH]) 11], allowing 
a smooth transition to the paraxial case, — > 0. 

It is worth noticing that the transformation to the he- 
licity basis is associated with the transition to the local 
coordinate frame with the z-axis attached to the current 
k-vector, which induces pure gauge Coriolis-type poten- 
tial A — —iWdkU, i.e., (A)jj = — *e* • (9k) Gj, where 
ei,2,3 = (e^,e^,K:) [13, [l^- At the same time, non- 
canonical operators and commutation relations ([2])- ([5]) 
essentially owe their origin to the projection onto the 
transversality subspace, which is equivalent to the diag- 
onalization of the potential A 10, 23]: Ab = dgA, i.e.. 



(Ab), 



-ie* • (^k) ejS, 



(9) 



While such diagonalization (which uncouples the two he- 
licity components) is an adiabatic approximation for a 
nearly-transverse paraxial wave beam propagating in an 
inhomogeneous medium [2^ [IJj , it is exact for transverse 
plane waves in free space where the helicities are truly in- 
dependent. 

The measurable expectation (mean) values of the 
0AM, SAM, coordinate, momentum, and energy ob- 
tained from the diagonal operators ((H])-® can be written 
as 



S = (E" (tk E"), 



L=(E' 



E^^-lE^ 



ctAr X k 



i9k 



E" 



E" 



crAf 



En, 

E" 



F ^ (E" k E"), W ^ (E" uj E' 



(10) 

(11) 
(12) 

(13) 
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Here = e^fc ^ (1 — cos 6') / aind, convolution implies 
summation over a — ±1 and integration in the k-space, 

and we assume normalization N ~ [e" E"'^ = 1 (see 

Appendix for details). While the SAM is purely intrin- 
sic (origin-independent), the 0AM, i n g eneral, has both 
intrinsic and extrinsic contributions '251: 



R X P 



(14) 



Equations ((T0)) - p4)) contain all the main observable 
results related to the AM and SOI of light. First, the 
CT-dependent non-paraxial Berry-phase term in L should 
be associated with the spin-to-orbit AM conversion 0, 
Hi, 12-14]. Particular cases of this term have appeared 
in [^,0,11. Second, the orbital and spin [iMll 

Hall effects of light are described by the two terms in the 
position of the center of gravity, Eq. ((T^ . Indeed, for 
a symmetric vortex beam propagating along the z-axis, 
the transverse coordinates of the center of gravity vanish, 
(A, y) = 0, after integration over 0, but any asymmetry 
of the field distribution along, say, the x-axis immediately 
causes an £- and cr-dependent shift along the orthogonal 
y-axis, Y together with tilt ^ (see the example 
in Section IV). 

We emphasize that our results (|6))-p4)) are exact and 
no approximations were made. They are equivalent 
to application of the canonical operators L, S, and f 
to the laboratory- frame field components {Ex, Ey, E^^)'^ 
supplied with the transversality condition. 



III. ENERGY FLOW APPROACH 



Remarkably, the same results, Eq. pl) )) - (|14l) . can be 
derived from an approach based on the separation of the 
spin and orbital parts in the Poynting energy fiow 0| . 
Let us consider a monochromatic beam-like field prop- 
agating in the positive z-direction. Pecularities of the 
(2 -|- 1)D formalism for such a problem are discussed in 
the Appendix, Eqs. (|I1)) - (IX9)) . 

The transverse center of gravity (jA9p obtained in the 
momentum representation from operator equation (|12|) 
can be equally derived from the traditional coordinate- 
representation definition 



R± (^) = ^J r±|E(ri,z)|'dV 



(15) 



where rj^ = {x,y), g — y^w/eo, and we used normal- 
ization / |E (r_L, z)!^ d^r_L = g"^ corresponding to A = 1 
(see Appendix). Substituting here Fourier representation 
()A6P with the helicity-basis expansion ()A2p . and using 
/r^e^'^^-^d^r^ = {27rf 6^ (kj.)5k, (k^ = (fc„fc^)) to- 
gether with expression ^ for the Berry connection, we 
arrive at Eq. (IA9|) which is equivalent to Eq. (IT2t . 



To derive the linear and angular momenta of the field, 
we use the Poynting vector which determines the momen- 
tum density (energy flow) p6l |: 



TT = 4 Im [E* X (V X E)] . 
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(16) 



Substituting here the Fourier decomposition (|A6p . we ob- 
tain 

TT = Re //e<k-k') r^.'*^.g.'* X X e-)^^, 
J J 27r 27r 

(17) 

where E"^' = E" (k'), e'^' = (k'), and summation over 
cr = ±1 is implied hereinafter. After some calculations 
the total momentum density PT)) can be decomposed into 
the orbital and spin parts as suggested in [g, tt = 

TT -|- TT : 



-^-■-E-E^'* (e^ • e'^'*) ^<L±^<L^ 
^ ^ 2 27r 27r 



(18) 

^ ^ 2 27r 27r ' 

(19) 

where = k ± k'. Note that in the decomposi- 

tion (|17p - p^ the two helicity components are exactly 
separated without an interference term [(|. Also, the 
"electro-magnetic democracy" discussed by Berry Q is 
accomodated, because switching to magnetic-field plane- 
wave helicity amplitudes, E'^ H'^ — ~iaE°', keeps 
Eqs. (IT7l)-([ni) invariant. 

The linear momentum of the beam (per unit z-length) 
is given by the 2D space integration of the momentum 
densities ([T7)) -(ITO1): P = / 7rd^rj_. In doing so, we find 
that the spin momentum density makes no contribution 
to the linear momentum [7]: P** = J tt^ (Pr± = 0, while 
the orbital contribution yields 



P° = P = / 7r°d^r_L = I k E-^ ^d^k± , 



(20) 



which obviously coincides with Eq. (|Tc 

The SAM and 0AM of the beam (per unit z-length) 
can be obtained by the 2D space integration of their den- 
sities, i.e.: 



L = / r X Tv°(frj_ , S = / r x n^'d^r^ 



(21) 



Substituting Eqs. (HH]) and into Eqs. ([HI) and em- 
ploying properties of Fourier integrals, we arrive at 



S = 













1 (e-xe-) 


E-' 




j OK 


E" 



(22) 
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L = / E''*e''* ■ {-ik X 9k) S'^e'^rf^kj. = / E"* {-ik x 9k) -B'"rf^k_L - / (t(Ab x k) E" 



d^k. 



(23) 



r 




FIG. 1: (Color online) (a) Bessel-beam distribution on 
the sphere in k-space with the azimuthal phase 2tt£. (b) 
Cylindrical caustic in the real space, an example of the closed 
orbit r on it, and the corresponding GO rays tangent to 
the caustic. Scalar phases are color-coded for £ = —1, 
do ~ 7r/4. Points A and A' on the caustic are connected by 
two paths: the straight line and the Poynting-flow helix. The 
phase matching yields the phase difference 2tt£ between the 
paths and quantization of the caustic radius. For circularly- 
polarized waves, the helical path brings about an additional 
Berry phase a$B, Eq. pH]! . 



where identity e*^* x e'" = ian and Eq. ^ were used. 
Clearly, the values of SAM and 0AM, Eqs. (HH) and dM]), 
derived from the Poynting energy flows are in perfect 
agreement with our operator formalism, Eqs. (1101) and 
([TT|) (see also Appendix). 



IV. APPLICATION TO BESSEL BEAMS 

Importantly, our theory has a number of directly ob- 
servable consequences. As the simplest example we 
take non-paraxial vector Besscl-bcam solutions which are 
eigenmodes of Jz constructed fro m p lane waves with well- 
defined helicity a (cf. [2, S HH, IIZI ) • The angular spec- 
trum of such beams is 



•Efi 



A" 5 {e - e^) e'^'l' , (24) 



where A"^ is a constant amplitude, 6q is the polar angle 
of conical distribution of the k- vectors, Fig. 1(a), and no 
summation over a is implied here. 

For the z-components of 0AM and SAM, Eqs. (HO]) 
and (ITT]) , or ([^ and (1^ . of a superposition of cr = ±1 
beams ([M]) we obtain ^281 ]: 



Herea = (|A+|^-|^-|^)/(|A+|^ 
aged hehcity and 



\A 



- |2n 



$B = J) As • dk = 27r (1 - cos ( 
c 



is the aver- 



(26) 



is the Berry phase associated with the contour C = 
{9 — 9o,(f> e (0, 27r)} formed by the k-vectors distri- 
bution on the sphere of directions, Fig. 1(a) 29]. The 
Berry phase is equal to the flux of the monopole field 
Fb — X Ab — k/k^ through the area of the k-space 
sphere bounded by the contour C. In this manner, the 
(T-dependent term in represents a monopole-flux con- 
tribution to the 0AM, cf. Eq. (87) in In the parax- 
ial limit the Berry-phase terms vanish as $3 ~ ttOq 0. 
The values (|25p evidence an apparent partial conversion 
from SAM to 0AM in non-paraxial light with the total 
AM being constant 0, Q , akin to the spin-to-orbit AM 
conversion upon focusing of polarized light M, M , Il2| - [l^ . 
Indeed, in the Richards- Wolf approximation [31| , the fo- 
cusing represents a geometric conical redirection of par- 
tial plane waves with their helicity being conserved. It 
is described exactly by the same transformation operator 
U{0, (f>) that describes transition to the helicity basis [3]. 

Simultaneously with a tr-dependent 0AM, the non- 
paraxial fields exhibit tr-dependent intensity distribu- 
tions related to the modified position operator. The real- 
space field of the circularly-polarized Bessel beam, calcu- 
lated via the Fourier transformation (IA6|) of Eq. 



(X A- I J, (0 + afee'(-+i)'^ J,+,+i (0 



where (p, ip, z) are the cyhndrical coordinates in real 
space, a — cos^ (^o/2), h = svc? (6*0/2), ^ = fc_Lp, fc_L = 
/c sin 6*0, fc|[ — fccos^o, Jn{£,) are the Bessel functions of 
the first kind, and the field components are written in the 

basis 

TT^' ^^v?^' ' '^'^^^ ^^^'^ ^ cylindrically 
symmetric intensity distribution, = given by 

cx (0 + (0 + 2a6j2^, (0] , (27) 

Below we show that the polarization-dependent intensity 
distributions ^} (see Fig. 2(a)) signify the SOI of light. 

The and cr-dependence of the radial intensity pro- 
file ([77)) can be explained via a geometrical-optics (GO) 
ray picture and the quantization of caustic underly- 
ing the maximum of the intensity. The rays associ- 
ated with a Bessel beam are those that form an an- 
gle 6q with the z-axis and touch a cylindrical caustic 
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6 



component = P^Yf = 
conversion is achieved at 6 



'jcotOoL^ [11]. The total 
0,7^1. 



V. CONCLUSION 

To summarize, we have revisited the problem of the 
identification of the spin and orbital angular momenta of 
nonparaxial light in free space. It has been shown that 
this issue is closely related to the determination of the po- 
sition of the center of gravity of a light beam or a wave 
packet. We have given an exact self-consistent solution to 
these problems in terms of quantum-operator formalism 
and using classical Poynting energy flows. In the helicity 
representation, taking into account the transverse nature 
of the electromagnetic fields, the operators of the 0AM, 
SAM, and position become diagonal, but exhibit non- 
canonical commutation relations. We have shown that 
the unusual features of these operators originate from 
the Berry-phase terms and can be associated with man- 
ifestations of the spin-orbit interaction of light. Indeed, 
anomalous Berry terms in the 0AM and position oper- 
ators describe spin-dependent part of 0AM (responsible 
for spin-to-orbital AM conversion) and spin-dependent 
shift of the center of gravity of light (i.e., the spin-Hall ef- 
fect of light). We have applied the general theory to sym- 
metric and asymmetric vector Bessel beams and found 
that our non-canonical operators indeed correspond to 
the observable quantities. The obtained Bessel-beam in- 
tensity distributions exhibit fine SOI splitting of caustics 
and Hall effects of light in perfect agreement with the 
derived 0AM and position operators. These effects can 
be observed experimentally in tightly focused fields. 

This work was supported by the European Commission 
(Marie Curie Action), Science Foundation Ireland (Grant 
No. 07/IN. 1/1906), the Australian Research Council 
(ARC), and von Humboldt foundation. We are grateful 
to A. Y. Bekshaev and M. V. Berry for fruitful discus- 
sions. 



Appendix A: Operator formalism for wave packets 
and beams 



One can separate two basic situations, for which the 
operator formalism of Section II can be adopted in a 
slightly different way. The first one is evolution of a wave- 
packet-like field localized in 3D space. Obviously, such 
field is nonmonochromatic and time-dependent. The 
plane-wave Fourier decomposition of the complex elec- 
tric field can be written as 



E 



(r, t) = J E (k) e*"— ('')*d3k , (Al) 



where d'^k — dk^dkydk^ = sin 6 dkdOdcf), luj(k) = fc is 
the dispersion relation, factor g — y^2uj/eo {eq is the vac- 
uum permittivity) is introduced for proper normalization 



of energy below, and the real wave electric field is given 
by £ (r, t) = Re E (r, t). In the helicity basis one has 



E (k) = E+ (k) e+ (k) + E- (k) e- (k) 



(A2) 



The energy of the wave-packet field is given by the 3D 
space integral of the intensity (we omit inessential con- 
stant factors) and can be written as: 



W = 



eo|Er + /io|Hnd^r 



E''*E''d^k= (E- 



E" 



(A3) 



Here H (r, t) = Re H (r, t) is the magnetic field, is 
the vacuum permeability, d'^r = dxdydz, summation over 
a — ±1 is implied hereinafter, and we performed some 
standard calculations with Maxwell equations and the 
Fourier transform (|Aip . Thus, the convolution implies 
3D integration of the field spectral amplitudes over the k- 
space. At the same time, to determine properly the state 

), one has to take into account the temporal 



vector 



dependence of the field, namely: 



E") =E'' (k) e 



— iC(;(k)t 



(A4) 



We assume normalization which has the meaning of 
the unit number of photons in the wave packet: = 

l^E'' E"^ = 1. Substituting the state vector (|X4)) with 

the definition of convolution (|A3I) into Eq. ([12]), we ob- 
tain the time-dependent position of the center of gravity 
of the wave packet moving in space: 



R (i) = - Im j E''*di,E''d^k - J aAg E" d^k + Vi, 

(A5) 

where the velocity of the wave-packet motion is given by 

V = y" (dwoj) E" ^d^k = j K E" ^d^k. 

Note that the same expression for the wave-packet cen- 
ter can be obtained by convolution of the canonical co- 

E'"\ = 



ordinate operator f = i(?k with the vector state 

E'^ (k) er^^i^)*- _ The Berry-connection term arises in this 
case from the 9k derivatives of the helicity basic vectors 
e"', Eq. The linear and angular momenta of the field, 
Eqs. pO| . (|lll) . and (IT3|) . can be calculated in a manner 
simflar to Eqs. (IX3)) - (|X5|) . 

The second typical problem that arises in optics deals 
with a beam-like monochromatic field (a; = fc = const) 
propagating in the positive z-direction and localized only 
in the transverse {x, y)-dimensions. In this case, it is 
natural to use the (2-1- 1)D version of quantum- like for- 
malism, where z instead of time plays the role of the 
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independent variable, whereas rj^ = (a;, y) is the ef- 
fective 2D space ahowing normalization of the trans- 
verse field distributions [3a]- Because of the monochro- 
maticity, only two components of the k-vector are in- 
dependent, and the z-componcnt can be expressed as 
kz — kz {'k±) = y^oj'^ — k'j_, k± = {kx,ky). This deter- 
mines the following 2D plane- wave Fourier decomposition 
of the complex time-independent electric field 36] : 

E (r_L, ^) = £: y E (k_L) e*k_L-r_L+»fc.(k^)^^2j^^ ^ ^^g^ 

where the real wave electric field is given by £ (r, t) — 
Re [E (r) e~*'^*] and the element of the 2D area of inte- 
gration is d^k_L = dkxdky. Alternatively, one can use 

E = E(0, ^) and (i^kj_ — k^ cos d sin Odd dtp in spherical 
coordinates with two independent dimensions (9, (p). The 
characteristic energy of the wave beam is, in fact, the en- 
ergy per unit z-length which is obtained by the 2D inte- 
gration of the time-averaged intesity over d^v±^ = dxdy: 



W=\j (£o]5]'+Mo]H]')rfV 



(A7) 



Here (r, t) = Re [H (r) e~*"*] , the overline stands 
for the time averaging, and we assumed the unit num- 
ber of photons per unit z-length in the beam: N = 



E" 



E" 



1. Thus, the convolution for beam-like 



fields implies 2D integration over (fc^, ky) or 



the 



k-space (these are equivalent unless we consider evanes- 
cent modes). To determine properly the state vector 

£'°'^, one has to take into account the z-dependence of 

the field, cf. Eq. (IX4l) : 



ik~ (k_|_)js 



(A8) 



E° 



Substituting definitions (|A7| and (|A8)) into Eq. ([12]), we 
obtain the z-dependent transverse position of the center 
of gravity of the propagating wave beam [36| : 

R_L (z) = -Im j E"* di^^E" d^\ii_ 

2 



(tAf 



E° 



d2kj_ + Vz 



(A9) 



where the 'velocity' of the motion along z is given by 



v = - / {d^^kz) 



E" 



d^k^ = / — 



E" 



d^k±. 



The linear and angular momenta (more precisely, their 
values per unit z-length) are calculated from Eqs. (jlOp . 
pT|) and in a similar manner (see also Sections III 
and IV). Note that despite the above (2 -|- 1)D quantum- 
like formalism, they are vectors in 3D space. This does 
not cause any difficulties if one uses z = idk^ = z which 
yields Z = z in. the 3D calculations, cf. [isj . 

It should be emphasized that, despite our using the 
same letters for the unifying formalism, quantities E, E, 
W, R, etc. have different meanings for the 3D-localized 
wave-packet polychromatic fields, Eqs. (|All) - (jA5p . and 
2D-localized monochromatic beams, Eqs. (jA6l) - (|A9p . 
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